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Abstract. We present a new chemical evolution model meant to be a first step in the self-consistent study of both optical 
and X-ray properties of elliptical galaxies. Detailed cooling and heating processes in the interstellar medium (ISM) are taken 
into account using a mono-phase one-zone treatment which allows a more reliable modelling of the galactic wind regime 
with respect to previous work. The model successfully reproduces simultaneously the mass-metallicity, colour-magnitude, the 
L x - L B and the L x - T relations, as well as the observed trend of the [Mg/Fe] ratio as a function of <r, by adopting the 
prescriptions of Pipino & Matteucci (2004) for the gas infall and star formation timescales. We found that a late secondary 
accretion of gas from the environment plays a fundamental role in driving the L x - L B and L x -T relations and can explain their 
large observational scatter. The iron discrepancy, namely the too high predicted iron abundance in X-ray haloes of ellipticals 
compared to observations, still persists. On the other hand, we predict [O/Fe] in the ISM which is in good agreement with the 
most recent observations. We suggest possible mechanisms acting on a galactic scale which may solve the iron discrepancy. 
In particular, mixing of gas driven by AGNs may preserve the gas mass (and thus the X-ray luminosity) while diluting the 
iron abundance. New predictions for the amounts of iron, oxygen and energy ejected into the intracluster medium (ICM) are 
presented and we conclude that Type la supernovae (SNe la) play a fundamental role in the ICM enrichment. SNe la activity, 
in fact, may power a galactic wind lasting for a considerable amount of the galactic lifetime, even in the case for which the 
efficiency of energy transfer into the ISM per SN la event is less than unity. 
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1. Introduction 

The most direct evidence for the presence of a non-negligible 
interstellar medium (ISM) in elliptical galaxies comes from X- 
ray observations. The presence of ~ 10 s - 1O IO M of gas in a 
hot (~ 10 7 K) X-ray emitting (L x ~ 10 39 - lO^ei-gs^ 1 ) phase, 
was first discovered by the Einstein satellite (e.g. Forman et 
' al. 1979) for ellipticals with a wide range of optical luminosi- 
ties. A correlation between X-ray (Lx) and optical (Lb) lumi- 
nosities is observed with roughly Lx oc L 2 B (e.g. O' Sullivan 
et al. 2001a), although the scatter is quite large. Such a trend 
can be explained by the fact that optically brighter galaxies 
are also more gravitationally bound (Ciotti et al. 1991). On the 
other hand, the X-ray emission from faint galaxies is domi- 
nated by discrete sources which have stellar origin (Brown & 
Bregman 2001) and the above correlation becomes approx- 
imately Lx oc Lb (Canizares et al. 1987; O'Sullivan et al. 
2001a). In what follows, we focus on massive systems which 
follow the Lx oc L 2 B trend and whose X-ray data tell us some- 
thing about their ISM properties. Concerning the scatter of the 
Lx - Lb relation, it is still not clear whether its origin can be 
traced to one of galactic environments (see e.g. the recent re- 
view by Mathews & Brighenti 2003 and references therein). 
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However, ellipticals located in denser regions do seem to ex- 
hibit higher Lx at a given Lb with respect to their counter- 
parts living in low density environments (Helsdon et al. 2001; 
Matsushita 2001, but see Matsumoto et al. 1997 and O'Sullivan 
etal. 2001a). 

As for clusters of galaxies, the X-ray properties of individ- 
ual ellipticals follow an Lx - Lx relation, although it is not 
a mere extrapolation from the relation observed on the inter- 
galactic scale (Matsushita et al. 2000), because the so-called 
'entropy floor' seems to be more active on a galactic scale. This 
fact can be easily understood in terms of the stellar feedback 
still ongoing due to Type la supernovae (SNe la) and low mass 
stars as well as other gas-dynamical processes (Matsumoto et 
al. 1997). A detailed study of these topics can be found in 
Kawata & Gibson (2003b), who found that radiative cooling is 
an important driver of this relation, although it does lead to late 
star formation (SF) and thus to galactic colours which are too 
blue with respect to observations. The additional energy feed- 
back provided by AGN may rectify this problem (e.g. Kawata 
& Gibson 2004). 

X-ray spectra also carry information pertaining to the 
chemical composition of the interstellar gas. The ASCA satel- 
lite provided the first reliable measure of the iron abundance in 
the hot ISM of ellipticals (e.g. Awaki et al. 1994; Matsumoto 
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et al. 1997). These observations led to the so-called iron dis- 
crepancy (Arimoto et al. 1997), the fact that the inferred iron 
abundance was much lower than the solar value, at odds not 
only with theoretical models for elliptical galaxies (Arimoto & 
Yoshii 1987; Matteucci & Tornambe 1987), which predicted 
that their ISM should exhibit [Fe/H]> 0, but also with the 
mean metallicity of the stellar component inferred from op- 
tical spectra. Arimoto et al. (1997) analysed different possi- 
ble sources of this discrepancy such as dilution by intracluster 
medium (ICM), different binary populations (and thus SNe la 
rates) with respect to spirals, the presence of dust, and the un- 
certainties involved in the modelling of the Fe L lines. In fact, 
Matsushita et al. (1997) showed that by taking into account sys- 
tematic errors in ASCA data, the iron abundance can possibly 
reach the solar value. This issue has been partly alleviated be- 
cause central iron abundances determined with ASCA did not 
take into account temperature gradients (Buote & Fabian 1998; 
Buote 1999), due to its limited spatial resolution. When a multi- 
temperature plasma model is adopted, the iron abundance in the 
hot ISM seems to increase from ~ 0.2 - 0.4 Fe at large radii, 
to solar (or slightly higher) inside a radius of ~ 50 kpc (see 
Mathews & Brighenti 2003). 

Moreover, thanks to XMM and Chandra satellites, it is now 
possible not only to obtain a more reliable measure of the 
iron abundance in the ISM, but also to observe other chem- 
ical species like O, Mg, Si, N, S, C (e.g Buote et al. 2003; 
Sakelliou et al. 2002; Gastaldello & Molendi 2002; Xu et al. 
2002; Matsushita et al. 2003). 

Elliptical galaxies have well-studied optical photo- 
chemical properties; it is natural to consider whether the the- 
oretical interpretations of these observable can simultaneously 
provide a satisfactory explanation for their X-ray properties. 
The X-ray spectra provide unique tools to study these reser- 
voirs of gas in which the evidence of the role played by re- 
cent SNe la activity and mass loss, due to stellar winds, can be 
found. 

Models for the chemical evolution of elliptical galaxies 
(e.g. Larson 1974; Matteucci & Tornambe 1987; Arimoto & 
Yoshii 1987) can successfully explain the mass-metallicity (e.g. 
Carollo et al. 1993; Gonzalez 1993; Davies et al. 1993; Trager 
et al. 1998, 2000a, 200ob) and the colour-magnitude (CMR: 
e.g. Bower et al. 1992a) relations, as well as the observed radial 
gradients in the line-strength indices and colours (e.g. Carollo 
et al. 1993; Peletier et al. 1990). Under the assumption that 
the SF efficiency increases with the galactic mass (Tinsley & 
Larson 1979; Matteucci 1994), and that the accretion timescale 
is lower in higher mass systems, Pipino & Matteucci (2004, 
PM04 hereafter) were able to match simultaneously the rela- 
tions presented above and the magnesium overabundance with 
respect to iron in the stellar component, which is observed to 
increase with galactic mass (e.g. Faber et al. 1992; Carollo et 
al. 1993; Davies et al. 1993; Worthey et al. 1992). 

The aim of this paper is to improve the detailed photo- 
chemical evolution code of PM04 with a new treatment of the 
hot ISM, in order to present, for the first time, a model able 
to predict in a self-consistent manner the properties of both 
stellar (i.e. optical) and hot gas (i.e. X-ray) components of the 
flux emitted by elliptical galaxies. We have modified the equa- 



tions of chemical evolution presented by PM04, by introducing 
a secondary gas inflow phase and explicitly taking into account 
the mass flow during the galactic wind. We next updated the en- 
ergetics by implementing the radiative cooling of the ISM fol- 
lowing Kawata & Gibson (2003a) and adopting the procedure 
described by Kawata & Gibson (2003b) to generate an X-ray 
spectrum of the model galaxies. In particular, we will show how 
the L x - L B relation depends critically on the interplay between 
the gas mass flow during the wind phase and the small sec- 
ondary accretion episode. The chemical code will be presented 
in Sec. 2, with particular emphasis on the novelties related to 
this new energetic formulation, as well as the method adopted 
in order to derive consistent X-ray properties. The models and 
the results will be shown and discussed in Sees. 3 - 6 and our 
conclusions will be summarized in Sec. 7. 

2. The model 

2.1. The chemical evolution code 

The adopted chemical evolution model is based on that pre- 
sented by PM04. In this particular case, however, we consider 
our model galaxies as a single zone extending out to 10/? e ff, 
where R e f[ is the effective radius (see Table 1), with instanta- 
neous mixing of gas. Moreover we take explicitly into account 
a possible mass flow due to the galactic wind and a possible 
secondary episode of gas accretion. Therefore, the the equation 
of chemical evolution for a single element i takes the following 
form: 



dGj(t) 
dt 



r M a m 

+ I tfr(t - T m )Q mi (t - T m )(p{m)dm 

nM B „ r ^0.5 

+A <f>(m)\ f(fi)Q mi (t - T mi )^j(t - T mi W 

r M B M 

+(1 - A) if/(t - T m )Q mi (t - T m )(f>(m)dm 

JM Bm 

+ I \fi{t - T m )Q mi (t - T m )(f>(m)dm 

Jm b„ 

dGdi) dGdi) 

+(— f^infall - WWi(t) + (— ^) acc , 

dt at 

where G,(f) is the normalized mass density of the element ; 
at the time t in the ISM. X,(f) is defined as the abundance by 
mass of the element i (see also Matteucci & Greggio, 1986, 
Gibson, 1997, and PM04 for a comprehensive discussion of 
this equation). The initial mass function (IMF) has the usual 
form (p(m) oc m~ (1+L35) , (Salpeter 1955), and it is normalized to 
unity in the mass interval 0.1 - 100M o . By means of this equa- 
tion we can calculate the evolution of 21 elemental species. 
The total ISM normalized gas mass density at each timestep 
is then simply G(t) = ll^G^t). The term -if/(t)Xi(t) gives the 
rate at which the element is subtracted from the ISM owing 
to the SF process (see Sec. 2.5 for details on the SF rate if/), 
whereas the integrals in the right-hand side of the equation give 
the rate at which the element i is restored into the interstellar 



dm 



(1) 



Pipino et al.: X-ray properties of massive elliptical galaxies 



3 



medium as unprocessed or newly-synthesized element by low- 
and intermediate-mass stars, and/or SNe la and Type II super- 
novae (SNe II). In fact, one of the fundamental points upon 
which our model is based, is the detailed calculation of the SN 
rates as well as the use of finite stellar lifetimes and detailed 
stellar nucleosynthesis (see Sec. 2.6). For type la SNe rate, we 
assume a progenitor model made of a C-O white dwarf plus a 
red giant and the formalism of Greggio & Renzini (1983) and 
Matteucci & Greggio (1986): 



R 



SNIa 



<KM B ) 

JM Bm Jfi„, 



(2) 



where Mb is the total mass of the binary system, Msm = 3M© 
and Mem = 16M are the minimum and maximum masses 
allowed for the adopted progenitor systems, respectively, [i = 
M2/MB is the mass fraction of the secondary, which is assumed 
to follow the distribution law: 



(3) 



Finally, fi m is its minimum value and y — 2 The predicted type 
la SN explosion rate is constrained to reproduce the present 
day observed value (Cappellaro et al., 1999), by fixing the pa- 
rameter A in eq. (2). In particular, A represents the fraction of 
binary systems in the IMF which are able to give rise to SN la 
explosions. In the following we adopt A = 0.09. 
On the other hand, the type II SN rate is: 

X16 
i/f(t - T m )(p(m)dm 



+ I \p{t - T m )<p(tn)dm , 
Ji6 



(4) 



where the first integral accounts for the single stars in the range 
8-16M , and Mjj is the upper mass limit in the IMF. The stellar 
lifetimes r„, are taken from Padovani & Matteucci (1993). The 
knowledge of the stellar lifetimes and the star formation history 
allow us to calculate the evolution of the stellar mass as given 
by: 



dM t (t) 
dt 



m(t) 



<K? - T m )(m - m rem )4>(m)dm , 



where m(f) is the smallest mass dying at the time t and m rem is 
the mass of the remnant of a star of initial mass m. Therefore 
we can easily evaluate the number of stars still alive, their dis- 
tribution in mass as well as their distribution as a function of 
the metallicity. In fact, for each simple stellar population (SSP) 
created, we store information on the chemical composition of 
the gas out of which this SSP has formed. 

The initial galactic infall phase enters the equation via the 
term: 



dGi(t) 

( ; )infall - ^i.infallCe 

dt 



(5) 



where Xi^u describes the chemical composition of the ac- 
creted gas, assumed to be primordial. C is the normalisation 
constant of the infall law, obtained by integrating this law over 
time and requiring that M\ um has been accreted at f gw (see 



PM04). After PM04, the initial collapse becomes faster as the 
galactic mass increases. 

At variance with previous semi-analytical works (PM04; 
Pipino et al. 2002, P02 hereafter, but see also Arimoto 1989; 
Gibson & Matteucci 1997), in which all the gas present at the 
time of the galactic wind is ejected into the IGM/ICM, in this 
model we take into account the possibility that only a fraction 
of the gas is able to escape the potential well. In order to do that 
we introduce the mass flow rate W(t) in the chemical evolution 
equation. The functional form of W(t) as well as its dependence 
on the energy content are given in Sec. 2.3. 

Finally, in order to explore the possible interaction between 
the pristine gas of the IGM/ICM and the hot ISM of the galac- 
tic models, we introduced an additional accretion episode, by 
means of the term (dGi(t)/dt) acc (see Sec. 2.4). 

The detailed photometric evolution for our model elliptical 
galaxies is obtained by applying the spectro-photometric code 
by Jimenez et al. (1998, see also PM04) to our SSPs. 

2.2. Energy evolution 

The chemical code presented in this paper features a new self- 
consistent energetic treatment which supersedes the previous 
one adopted by PM04. In particular, calling u t h the thermal en- 
ergy per unit mass, the equation governing the energy evolution 
of the gas is: 



diixh 
dt 



- SN 



M, 



n z H A(Z,T) 1 d M 

+ — — ■ («vir - «th) , (6) 



M gas dt 



where the first term on the right-hand side is the heating rate 
per unit mass (Sec. 2.2.2), while the second is the cooling rate 
per unit mass (Sec. 2.2.1). The last term represents the varia- 
tion in the energy content due to the stellar mass loss (injected 
into the ISM and being at the virial temperature of the system, 
see eq.|5}. This is a clearly simplified case with respect to those 
presented in other works (e.g. Ciotti et al. 1991, Loewenstein 
& Matthews, 1987, Hattori et al. 1987). Nevertheless, because 
we cannot treat gas dynamics, we believe that eq. (J6J is a good 
starting point for a detailed inclusion of stellar feedback in a 
chemical evolution model. The temperature of the ISM is re- 
lated to the energy per unit mass by means of: 



«th - 



-kT 



Equation (jS) has the form: 



du\<a 
~dt~ 



= /(«th) ■ 



(7) 



(8) 



Therefore, in order to solve equation ijfji for each time-step, we 
employ an iterative process. In particular, we implemented a 
routine based on Brent's Method (see Press et al. 1986) in the 
chemical code. 

This treatment is needed in order to study the evolution of 
the thermal budget of the ISM itself. The previous formalism, 
based on the thermal energy of the supernova remnant (SNR), 
is now part of the new procedure, so that we can distinguish be- 
tween the energy content of the gas and the energy trapped in 
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the hot bubble of the SNR. The latter is eventually restored to 
the ISM when the SNR merges with the surrounding medium. 
Moreover, the inclusion of a consistent calculation of the cool- 
ing functions allows us to compute the X-ray flux emitted by 
the model galaxies (see Sec. 2.7). 

Equation (J6J relies on the assumption that each gas parti- 
cle shares the same energy, whether it is gas being used to form 
stars or gas escaping in the galactic wind. This limitation is due 
to the single-phase nature adopted for the ISM. A more sophis- 
ticated treatment would require a two-phase model in which 
mass exchange is allowed by evaporation (cold phase— >hot 
phase) and cooling (hot— >cold) (Harfst, Theis & Hensler 2003). 
The model presented in this paper, however, is the first step 
in modeling both the X-ray and optical properties in a self- 
consistent manner. 

The relation: 

kT vir = ^^, (9) 

l « v ir 

(e.g. Romano et al. 2002) gives us the initial value for the en- 
ergy per unit mass; R vn - is the virial radius for a spherical sys- 
tem of mass M v ; r , namely the mass of luminous and dark matter 
inside 10 Tiff- 
in the case of secondary accretion, we take into account that 
the infalling gas has a mean energy assumed to be equal to the 
virial temperature (according to eq. |9j, therefore influencing 
the final ISM temperature. This can be easily done by adding 
the term -^p • (u vlr - u±) in the right-hand side of eq. (6). 

We limit our calculation to temperatures above T = 10 4 
K, since the cooling functions become more complex at lower 
temperatures (e.g. Omukai 2000). This condition, however, did 
not affect the final results, since in all the runs completed the 
gas temperature was always higher than this minimum value. 

2.2.1 . The cooling term 

After Kawata & Gibson (2003a, where we address the reader 
for a more detailed description), we evaluate the cooling rate 
A(Z, T) as a function of gas density and metallicity by means 
of the cooling curves computed by MAPPINGS III (Sutherland 
& Dopita 1993). In particular, we used a bi-linear interpolation 
in [Fe/H] and T in the range [-3,0.5] in metallicity and [10 4 K, 
10 9 K] in temperature. Since our models span a wider range 
in [Fe/H] during their evolution, we assumed that the cooling 
rate of the gas with metallicity [Fe/H]< -3 was the same as 
for gas with [Fe/H]= -3 in order to avoid uncertainties due to 
extrapolation. A similar assumption has been made for gas with 
metallicity [Fe/H]> 0.5. 

2.2.2. The heating term 

The term describing the heating of the gas in equation is 
defined as: 

FsN = eSNIa^O^SNIa + fSNII^O^SNII > (10) 

where Eq = 10 51 erg is the assumed initial energy budget for the 
SNR, esNia (esNii) and tfsNia (^snii) are the efficiency of energy 



transfer and the explosion rate for SNe la (SNe II), respectively 
(see eqs. 2 and 4). 

At variance with previous works (P02; PM04), here we 
adopt a parametric value for both esNia and esNii- This allows 
a more flexible treatment of the heating term in order to test 
different scenarios. In particular, before the galactic wind, we 
focus on three different cases, namely esNia = fSNii =0.01, 
esNia = esNii =0.1 and e S Nia = csnii =1. We note in passing 
that our fiducial case (esNia = esNii =0.1) is in good agreement 
with the results of Thornton et al. (1998), who modelled SNR 
evolution within a range of different environments finding that 
the fiducial value adopted here can be considered as typical. In 
the following, the fiducial case has been adopted unless oth- 
erwise noted. Furthermore, we have done preliminary calcula- 
tions for typical values of metallicity, density and temperature 
relevant to our models, based on the Cioffi et al. (1988) for- 
malism (see also Gibson 1997). In particular, considering that 
a SNR merges with the ISM (and hence eventually shares its 
energy with the surrounding medium) when its expansion ve- 
locity equals the gas sound speed, we found that the typical 
fraction of the initial energy budget available for the ISM is 
only a few percent, again in agreement with our fiducial case. 

After the galactic wind has developed, the gas content of 
the galaxy decreases quite sharply, whereas the temperature is 
almost constant or increases. When we take into account these 
new conditions in the Cioffi et al. (1988) equations, we have 
an increase in esNia up to values as high as 0.5-0.7 (see Table 
1). This is the most natural way to inject energy and drive a 
wind in our model, but, in a more general treatment (part of) 
this energy could be also provided by an AGN. 

A more detailed and self-consistent treatment of the SNR 
evolution as a function of gas metallicity and temperature will 
be included in a forthcoming paper. A more general form for 
the heating term would request the inclusion of the energy input 
from stellar winds (see e.g. Loewenstein & Mathews, 1987). Its 
typical value for massive stars is ~ 1 % of the assumed SNII en- 
ergy budget (Mathews & Baker, 1971). Therefore the inclusion 
of stellar winds as a heating source has a negligible impact on 
the final results (Bradamante et al. 1998; Gibson 1996). 

2.2.3. The development of the galactic wind 

We define the gas binding energy E^m (see PM04 for its ana- 
lytic definition) as the work required to carry a gas particle out 
to 10^ e ff. According to PM04, we assume that the dark mat- 
ter (DM) is distributed in a diffuse halo ten times more mas- 
sive than the baryonic component of the galaxy with a scale 
radius Rd- d ± = 10R e ff- The DM profile is taken from Bertin et 
al. (1992). 

In this paper we define the time when the galactic wind 
onsets (f gw ) as the solution of the following equation: 

£th(fgw) = £Bin(fgw) , (11) 

here the thermal energy of the gas is simply described by 

E th = MthM gas . (12) 

For every time-step in which the thermal energy is greater than 
the binding energy, we set the energy loss in the wind as the 



Pipino et al.: X-ray properties of massive elliptical galaxies 



5 



value given by AE = E± - E^m, whereas in previous works 
(P02; PM04) all the thermal energy was allowed to escape. 

Before introducing the modelling of the SF history and the 
possible exchanges of matter with the surrounding medium at 
t > f gw , we spend a few words on the differences between this 
new treatment for the energy evolution and our previous works. 
Since the heating term is not dramatically different from the 
energy source powering the winds in P02 and PM04 models, 
which featured an average SN efficiency of ~ 20%, the most 
relevant difference is the cooling of the ISM. This results in a 
lower effective contribution by a single SN to the energetic bud- 
get and, thus, delayed galactic winds with respect to previous 
models. Moreover, owing to the one-zone nature of this model, 
the adopted galactic radius R becomes an important parame- 
ter in determining the cooling term, because the latter scales 
as Pg as (and thus as R 6 ). The influence of the adopted radius 
on the ISM X-ray properties will be discussed in Sec. 5.2. The 
detailed energetic treatment will provide us with the tools to 
evaluate the amount of gas which has enough energy to escape 
from the galactic potential well and to control whether the con- 
ditions for the SF hold during the different stages of the galactic 
evolution. 



2.3. The mass flow rate 

The mass flow rate in the right-hand side of equation is 
defined as: 

AM 1 



W(f) = 



(13) 



At M lum 

Given the assumptions leading to eq. (6) and the fact that the 
amount of energy which can drive the wind is AE, the mass 
of gas which can blow away from the galactic potential well at 
each time-step is given by: 
AM _ AE 
E t h 



(14) 



gas 



Therefore, the energy feedback enters in the equation of the 
chemical evolution by means of eq. (13). 

2.4. Second Infall Episode 

Theoretical models based on the Cold Dark Matter paradigm 
predict that massive halos keep accreting either dark or bary- 
onic matter all their life (e.g. Voit et al. 2003). Since our aim is 
to reproduce the properties of very bright ellipticals, we further 
modified our chemical evolution code by adding the possibility 
of a late time accretion of gas. In this way we partly relaxed the 
somewhat rough approximation that the infall suddenly stops 
at the occurrence of the galactic wind. It is likely, in fact, that 
even in the wind phase, small amounts of gas can flow inward, 
e.g. in the presence of an inhomogeneous ISM. On the other 
hand, the optical properties of elliptical galaxies require that 
the main accretion episode must have been very fast compared 
to the galactic lifetime and that most of the galaxy has to be 
assembled during this stage (e.g. PM04; Ferreras & Silk 2003). 
Therefore, in the following we shall refer only to the main ac- 
creting episode as the infall. In order to handle easily the pa- 
rameters governing the secondary inflow, we chose to model 



independently these two accretion episodes and we adopted the 
following expression for the second one: 

M acc 



dGi(t) 

( T )acc — 

at 



(15) 

namely we assume that the model galaxy accretes a gas mass 
M acc in a uniform way all along its lifetime (to). This choice can 
be motivated by the fact that the secondary accretion is com- 
pletely different from the fast initial infall. Eq. (15), in fact, 
models in a simple way interactions with the environment or 
random accretions of a singular gas cloud occurring on long 
timescales, within a one-zone chemical model with instanta- 
neous mixing. Also, in this case the accreted gas is assumed 
to be primordial and its temperature is given by equation l|9}- 
In Sec. 5 we show the importance of this secondary accretion 
episode, in governing the late time behaviour of the galaxies. 



2.5. Star formation history 

The SF rate (SFR) is assumed to be proportional to the gas 
density, according to: 



iff(t) = v • p gas (0/Pgas,0 - 



(16) 



The SF efficiency v is taken to be an increasing function of 
the galactic mass following the prescriptions of the PM04 best 
model. 

An improvement of this model with respect to PM04 is that 
we can keep the gas temperature and its cooling time under 
control. These quantities allow us to check in a self-consistent 
manner whether the SF process might continue or not. In par- 
ticular, the SF takes place if the condition t coo \ < ts is satisfied 
by the ISM gas. The cooling timescale is defined as: 

Eth 3 Pgas kT 



^cool 



Sjit 2fim p n 2 H A(Z, T) 



While the free-fall time is: 

3?T 



t S = ( 



32Gp g as 



(17) 



(18) 



As it can be seen from Fig. 1, the condition is satisfied at 
the beginning of the SF process. Since then, because ellipti- 
cal galaxies cannot stop or regulate the SF process (Elmegreen 
1999), the model galaxies undergo a rapid and very efficient 
conversion of gas into stars until the galactic wind occurs. On 
the other hand, after the galactic wind the gas is so hot and rar- 
efied that the f coo i < ts condition cannot hold anymore, there- 
fore leading to a sudden stop in the SF, in agreement with pre- 
vious works (PM04). 

Finally, when the galactic wind stops, late SF episodes are 
still inhibited by the high temperature and very low density of 
the ISM, which imply f coo i > ts- It is worth noting, however, 
that the chemical evolution equation presented here holds even 
in the case in which wind and SF are present at the same time. 

2.6. Stellar Yields 

In order to follow the evolution of the chemical elements, we 
utilize the following nucleosynthesis prescriptions. 
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Fig. 1. Star formation rate (upper panel), f coo i, tg (lower panel) 
versus time for model Lai (see text and Table 1 for details). 

1. Yields of van den Hoek & Groenewegen (1997, as a func- 
tion of metallicity) for single low and intermediate mass 
stars (0.8 < M/M Q < 8). 

2. Yields by Nomoto et al. (1997) for SNe la. 

3. Yields of Thielemann et al. (1996), which refer to the solar 
chemical composition, for massive stars (M > 8M Q ). 

2. 7. X-ray properties 

In order to study the X-ray emission from model galaxies in a 
detailed manner, taking into account the line emissions and the 
energy range of the actual observed spectra, we applied the pro- 
cedure described by Kawata & Gibson (2003b) to our chemical 
code, which we briefly outline here. In particular, using the el- 
emental abundances in the ISM and the gas temperature and 
density as provided by the chemical code, we create a spec- 
trum in the energy range 0.1-20 keV, with 1800 bins, by means 
of the vmekal plasma model (Mewe et al. 1985; Mewe et al. 
1986; Kaastra 1992; Liedahl et al. 1995). Then we generate a 
fake observed spectrum by convolving the spectrum with the 
epn_ff20_sY9_thin.rmf response function for the EPN detector 
on-board the XMM-Newton satellite in XSPEC ver. 11.1.0 en- 
vironment. The exposure time is 40 ks and neither absorption 
component nor background were taken into account. The spec- 
trum is rebinned in order to have at least 25 counts per bin, and, 
finally, fitted with XSPEC vmekal model in the same energy 
range of the different observations. We stress that the abun- 
dance obtained as free parameters in the fitting procedure, are 
in good agreement with the output of the chemical code. 

A detailed analysis of this procedure, as well as a compari- 
son with other works in which the X-ray emission is modelled 
as simply proportional to p 2 T 1/2 , can be found in Kawata & 
Gibson (2003b). At variance with that paper, however, because 



we are presenting a one-zone model, we will not take into ac- 
count aperture effects. 

3. Galactic models 

We run models for elliptical galaxies with baryonic mass 
M\ um = 10 11 (models labelled as L) and M\ um = 1O 12 M (mod- 
els labelled as H). M\ um is the 'nominal' mass of the object, i.e. 
the mass of gas accreted and possibly turned into stars during 
its active evolution (recall that we stop the infall when t ~ f gw ). 
Furthermore, we name the models '0', 'al' or 'alO' if they have 
no secondary accretion, accrete 1 per cent or 10 per cent of 
Mi um , i.e. M acc =0, 0.01Mi um and 0.1Mi um , during their whole 
evolution, respectively. Models called HSN1 and HSN100 are 
similar to model Hal, except for the SN efficiency (1 and 100 
per cent, respectively). Since the secondary accretion episode 
implies a very mild accretion rate, the SF history, and hence the 
photochemical properties are not affected by it; only the late 
behaviour of the ISM properties depends strongly on the sec- 
ondary accretion, therefore affecting the X-ray luminosity and 
temperature. A summary of the input parameters are shown in 
Table 1, where model name (Col. 1), mass (Col. 2), effective 
radius (Col. 3), SF efficiency (Col. 4), infall timescale (Col. 
5), SNe efficiency (Col. 6-8), accreted mass in the secondary 
episode (Col. 10) as well as the times at which the galactic 
wind starts and stops (Col. 9 and 11, respectively) are given. 

4. The stellar component - photochemical 
properties of the models 

As it can be seen from Table 1, in this paper we follow the 
prescriptions given by PM04 in order to reproduce the major- 
ity of photochemical observables for elliptical galaxies, namely 
that the SF efficiency has to increase with the galactic mass, 
whereas the infall timescale must decrease with luminosity. It 
is worth noting that the values for v and r are very close to those 
presented by PM04. Therefore we expect that our model suc- 
ceeds in reproducing the mass-metallicity and the [Mg/Fe] -cr 
relations as well as the CMRs. We recall here that the observa- 
tions of the increase of line-strength indices (such as Mg2 and 
(Fe)) with galactic velocity dispersion (e.g. Carollo et al. 1993; 
Gonzalez 1993; Davies et al. 1993; Trager et al. 1998, 2000) 
is generally interpreted as a metallicity sequence in which the 
stars of more massive galaxies are richer in metals with re- 
spect to the stellar component of less massive objects. This 
metallicity sequence also explains the slope and tightness of 
the observed CMR (e.g. Kodama & Arimoto 1997). On the 
other hand, the observed increase in the mean stellar [Mg/Fe] 
in the central region for ellipticals with galactic velocity dis- 
persion (e.g. Faber et al. 1992; Carollo et al. 1993; Davies et 
al. 1993; Worthey et al. 1992, but see Proctor et al. 2004), 
tells us that the most massive galaxies must have formed in 
a shorter timescale compared to the less massive ones (time- 
delay model, see Matteucci 2001). We address the reader inter- 
ested in a detailed discussion on the topic and in the role played 
by v and r in reproducing the above relations to Matteucci 
(1994) and PM04 (see also Thomas et al. 2002; Ferreras & Silk 
2003). 
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Model 




«eff 


V 


T 


e SNII 


f SNIa 


f SNIa 






wind stop 




(Me) 


(kpc) 


(Gyr 1 ) 


(Gyr) 




« < t„) 


(t > ; pw ) 


(Gyr) 




(Gyr) 


LO 


10" 


3 


15 


0.3 


0.10 


0.10 


0.50 


0.82 





never 


Lai 


10 11 


3 


15 


0.3 


0.10 


0.10 


0.50 


0.82 


0.01 


7.8 


LalO 


10" 


3 


15 


0.3 


0.10 


0.10 


0.50 


0.82 


0.1 


3.9 


HO 


10 12 


10 


25 


0.2 


0.10 


0.10 


0.70 


0.49 





5.2 


Hal 


10 12 


10 


25 


0.2 


0.10 


0.10 


0.70 


0.49 


0.01 


3.7 


HalO 


10 12 


10 


25 


0.2 


0.10 


0.10 


0.70 


0.49 


0.1 


2.6 


HSN1 


10 12 


10 


25 


0.2 


0.01 


0.01 




no wind 


0.1 




HSN100 


10 12 


10 


25 


0.2 


1.00 


1.00 


1.00 


0.13 


0.1 


never 



Our model predictions for the photochemical properties re- 
lated to the optical part of the spectrum are summarized in 
Table 2 and compared with observation in Fig. 2 and Fig. 3. 
In this section, however, we simply refer to model L (which 
represents models L0, Lai, LalO) and model H (which rep- 
resents models HO, Hal, HalO), since the optical properties 
are not influenced by M acc . In particular, in Table 2 we show 
the predicted line-strength indices (see below) in Cols. 2-5, the 
[(Mg/Fe)] in Col. 6 and the luminosities of the whole galaxies 
in the K and L-band in Cols. 7 and 8, respectively. The pre- 
dicted indices are evaluated by calculating the mean [Mg/H], 
[Fe/H] and [Mg/Fe] in the stellar component of the galaxies, 
and then converted into Mg2 and (Fe) by means of a suitable 
calibration. In particular, the indices obtained by Matteucci et 
al. (1998) with the Tantalo et al. (1998) a-enhanced SSPs are 
labelled with T, whereas W refers to the Worthey (1994) SSPs 
(see PM04). For a detailed discussion on the adopted proce- 
dure and a comparison between the two sets of calibration see 
PM04. 

Before discussing the agreement with the observed mass- 
metallicity relation, we stress that the observed ones refer to the 
central (i.e. ~ 0.1/? e ff) galactic region, whereas we are present- 
ing the results for a one-zone model extending out to 10/? e ff- 
We also recall that the strength of the metallic lines is different 
on these two scales, being weaker at larger radii (e.g. Carollo et 
al. 1993). Therefore, in order to compare our predictions with 
observations in Fig. 2, we must convert our average results for 
the whole galaxy (entries of Table 2) into a precise estimate for 
its central region. In particular, we derived the line-strength in- 
dices for the central zone as Indexo = Index{\R e s) + AIndex, 
where AIndex is the absolute value of the theoretical gradients 
predicted by PM04 best model within 1 . In order to do that, 
we assume that Index(\R e s) ~ Index(lOR e g), namely that the 
mean stellar metallicity, either measured or predicted, at ~ R e g 
is representative of the whole galaxy (e.g. Arimoto et al. 1997). 
It is worth noticing that when considering line-strength indices 
derived by a-enhanced tracks (solid curves in Fig. 2), the nega- 
tive gradient in [Fe/H] is balanced by a positive radial gradient 
in [Mg/Fe] (see PM04), so the values for Mg2 and (Fe) pre- 
dicted by the one-zone model not only preserve the slope, but 
also lie very close to the average observed ones. 

Our fiducial models can reproduce reasonably well also the 
[Mg/Fe] -cr relation, once the entries in Table 2 are corrected 
for the radial gradient. In particular, we assume the values 
given by PM04 for their case lib, namely A [Mg/Fe] =0.326 for 
model L and A [Mg/Fe] =0. 169 for model H, leading to the pre- 



v 



c: 



cio o 

s 




H h 



H 1 1 h 




10 11 12 

Log (M tot /M ) 

Fig. 2. Predicted mass-metallicity relation for the central zones 
of our models. The curves are obtained by correcting the aver- 
age properties of the whole galaxy with the theoretical metal- 
licity gradients by PM04 (see text for details). Results for a- 
enhanced calibration (solid line) and not a-enhanced calibra- 
tion (dashed line) are shown versus a collection of data from 
Kuntschner (2000) and Kuntschner et al. (2001). 



dieted central values [Mg/Fe] =0.1 80 (L) and [Mg/Fe] =0.392 
(H), which are in good agreement with the observations. 

We show the predicted V - K versus M v and U - V versus 
M v relations in Fig. 3 for two different ages for the galactic 
models (the fiducial case with esN = 0.1) compared with the 
data of Bower et al. (2002b). In order to satisfy the condition 
given by the CMRs, we assume that the galaxy in model H 
forms at redshift 3 (age of 1 1.3 Gyr) while in model L it forms 
at redshift 2.6 (age of 1 1 Gyr) in an Q m = 0.3, Q A = 0.7, 
H = 70kms _1 Mpc 1 cosmology. The fact that the less mas- 
sive galaxies might be also the youngest, has been recently sug- 
gested by several authors from the analysis of their observed 
line-strength indices (e.g. Thomas et al. 2002). 

In order to show how the choice of the age can affect the 
predicted CMRs (age-metallicity degeneracy, e.g. O'Connel 
1976), we show a case in which the galaxies in both models 
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Table 2. Predicted optical properties of the model galaxies 



Model 


Mg 2 




<Fe> 




KMg/Fe)] 


Lb 


U 




T 


W 


T 


W 




(10 10 L o ) 


(10 10 L o ) 


L 


0.301 


0.232 


3.00 


2.68 


0.506 


1.2 


7.7 


H 


0.304 


0.220 


2.93 


2.57 


0.561 


11.2 


80 


HSN1 


0.431 


0.350 


3.94 


3.78 


0.462 


18 


116 


HSN100 


0.268 


0.172 


2.61 


2.10 


0.689 


4.4 


19.8 



The indices obtained with the Tantalo et al. (1998) SSP are labelled with T, whereas W refers to the Worthey (1994) SSP 
(see text). The values presented here refer to the whole galaxy. 



> 



> 




5. The gaseous halo 
models 



X-ray properties of the 



Fig. 3. Predicted CMRs for our models and different redshift 
of formation compared with the fit to the Bower et al. (2002b) 
data (solid line). A measure of the observational scatter is given 
by dashed lines. 



L and H form at redshift 5 (i.e. they have an age of 12.3 Gyr 
for the assumed cosmology), with full triangles in Fig. 3. 

Either models with a very low (HSN1) or maximum 
(HSN100) SN efficiency fail in reproducing the above relations 
(Figs. 1 and 2). As expected, model HSN1 hardly develops a 
wind, thus allowing SF to take place until the present time. SNe 
II explosion rate inconsistent with observations (Cappellaro et 
al. 1999), too blue colours and, at the same time, a too high 
stellar metallicity with respect to model H (and hence to the 
observed values) are natural consequences of this particular SF 
history. 

Model HSN100, on the contrary, undergoes a galactic wind 
too early, thus not allowing for the necessary enrichment of 
the stellar component up to (Z) ~ 0.02. As it can be clearly 
seen from the entries of Table 2, in this case the line-strength 
indices are too weak, and the galaxies exhibit very low lumi- 
nosities. For these reasons, model HSN1 and HSN100 will not 
be analysed in further detail in the remainder of the paper. 



5.1. The energy evolution of the ISM 

Once we have tested our models in their photochemical proper- 
ties, showing that the new energetic formulation preserves the 
results achieved by PM04, we dedicate this section and the fol- 
lowing one to our new results. In particular, we first concentrate 
on the temporal behaviour of cooling, gas density and temper- 
ature of the ISM, in order to understand what determines the 
L x - L B and L x - T x relations. Then we shall analyse the chem- 
ical composition of the predicted X-ray haloes. A summary of 
the predicted properties for the ISM of the model galaxies is 
shown in Table 3. For each model we shall show the role played 
by secondary accretion in determining the final hot gas mass 
(Col. 2), chemical abundances (Cols. 3 and 4), X-ray luminos- 
ity (Col. 5) and, finally, temperature (Col. 6). 

As it can be clearly seen from the entries of Table 3, mod- 
els which accrete ~ 0.1Mi um or more during their entire life- 
time, undergo strong cooling and exhibit too low temperatures 
to emit in the X-ray. In any case, a variation in M acc of a few 
percent can dramatically change the late time evolution of the 
galaxies. This finding can be explained by means of Fig. 4, us- 
ing model H as an example. First of all, it is worth noticing 
that the three models show no relevant differences in the heat- 
ing source, because at late times only SNe la are exploding and 
their rate is a function of the SF history (which is the same, 
since the secondary accretion affects neither the SFR nor the 
galactic wind occurrence). Model HO undergoes galactic wind 
for a longer period. At every time-step the condition > E^m 
is satisfied and this implies AE > and AM > 0, which, in turn, 
leads to an increase in the heating per unit mass and a decrease 
in the cooling (due to the decreasing density). This mechanism 
leads to slightly high temperatures (> 2 keV) and devoids the 
galaxy of most of its gas. The opposite behaviour is showed by 
model Ha 10, where the condition E± > Esm cannot hold for a 
long time interval, owing to the higher mass accretion rate and 
faster cooling with respect to model HO. It should be remarked 
that, in models Ha 10 and La 10, the accretion rate at t > t gw ex- 
ceeds the mass loss rate by a factor of 2, although M acc is much 
lower than dM w /dt integrated over the entire galactic lifetime. 
This is due to the fact that, ~ 100 Myr after t gw , only low mass 
stars (i.e. mass < 3M ) are still alive in the galaxy, and their 
contribution to the returned mass is quite low. When the galac- 
tic wind stops, all the gas either ejected by stars or accreted 
from the ICM is retained and leads to an even stronger cooling. 
The galaxy eventually exhibits a cold ISM. 
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Table 3. Predicted X-ray properties of the model galaxies 



Model 


Af gas /Mi um 


[Fe/H] 


[O/Fe] 


Lx 


kT x 








/ 1 A41 ,,-1 \ 

( l u erg s ) 


(keV) 


LO 


0.14- 10~ J 


0.85 


-0.56 


9.3 ■ 10- 3 


1.2 


Lai 


0.014 


0.85 


-0.68 


2 


0.5 


LalO 


0.15 


0.92 


-1.05 




0.001 


HO 


o.3 ■ ur 2 


1.23 


-0.87 


16.7 


1.7 


Hal 


0.05 


1.3 


-1.04 


183 


1 


HalO 


0.17 


1.60 


-1.17 




0.001 


Hal-0.3Zo 


0.05 


1.4 


-1.02 


175 


0.98 


Hal-Z 


0.05 


1.3 


-0.98 


254 


0.84 



Models Hal-0.3Z o and Hal-Z G are similar to model Hal, but the second accretion episode has a metallicity of 0.3Z e and 
Z e , respectively. 
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Fig. 4. Temporal behaviour of the ISM iron abundance, tem- 
perature, gas density and cooling term predicted for models HO 
(solid lines), Hal (dotted) and HalO (dashed). 



Models with M acc ~ 0.01M mm (i.e. Lai and Hal), on the 
contrary, show a reasonable behaviour of the gas temperature, 
with a good balance between heating sources and cooling term, 
given the fact that the accreted gas has the virial temperature of 
the galaxy. Similar conclusions were reached in previous works 
(e.g. Ciotti et al. 1991; Lowenstein & Mathews 1987; Ferreras 
et al. 2002). Ciotti et al. (1991), pointed out the importance 
of a late time inflow phase in order to explain the high X-ray 
luminosities of the most massive objects, although the gas ac- 
cretion they considered was limited only to the innermost re- 
gions of the galaxies and the inflowing gas had a stellar origin. 
From these results, in the following we focus on models Lai 
and Hal. 

In Fig. 5 we show the spectrum predicted for model 
Hal following the procedure described in Sec. 2.7. It is in 
good qualitative agreement with the observed spectra (see e.g. 
Gastaldello & Molendi 2002), and we notice a strong feature, 
centered at 1 keV, due to the very high Fe abundance predicted 
by our models. Among the metal lines, in fact, the Fe L shell 
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Fig. 5. Predicted X-ray spectrum for model Hal. For the 
adopted procedure see Sec. 2.7. 

gives the highest contribution, owing to the fact that the ISM of 
the model galaxies is rich in iron. In order to have more quan- 
titative constraints for our model (with particular emphasis on 
models Lai and Hal), we extract information about L x , T x and 
[Fe/H] from the fake spectra and compare them with observa- 
tions. 

5.2. The L B - L x relation 

First of all, we show how our models can match the L B - L x 
relation, which is the best determined link between the opti- 
cal and the X-ray part of the spectrum (either theoretically or 
observationally, see Sec. 1). In Fig. 6 we show our predictions 
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for model Lai and Hal (filled triangles), compared with the 
mean relation derived from the large sample of O' Sullivan et al. 
(2001a, solid line). Our points lie well within the observed re- 
gion (a measure of a 1 cr deviation is given by the dotted lines), 
although the predicted slope seems to be slightly flatter with 
respect to the fitted relations. 

An interesting result is that the scatter in the observed 
L B - L x relation can be totally explained either by small dif- 
ferences in the accretion history or variations in density (as we 
explain later in this section, see also Ciotti et al. 1991). In order 
to show this, we plot in Fig. 6 models with M acc = 0.005Mi um 
and M acc = 0.007Mi um (filled squares and stars, respectively) 
as well as models L0 and HO (empty squares). The difference 
in the predicted L x at a given L B for the accretion models and 
the cases without secondary inflow spans the same region that 
is between the fit to the observed data (solid line in Fig. 8) 
and the lcr boundary (dotted line) in the case of the high mass 
model. For the low mass one, instead, the range in the pre- 
dicted values covers a much larger region. This fact is a con- 
sequence of the adopted formalism for the energetics. In fact, 
with the exception of the very beginning of the galactic wind, 
when £ t h » £Bin, the condition on the mass flow presented in 
eq. (14) coupled with a mild accretion of slightly colder exter- 
nal gas, implies that E± > E^m at each time step. In this cases, 
the mass accretion rate is 0.2 - 0.5 times the value of the mass 
loss rate. Therefore it gives a non-negligible contribution to the 
building-up of the X-Ray emitting ISM. Even small changes in 
the accretion rate can cause differences in the predicted proper- 
ties. The higher density of the low mass model with respect to 
model H amplifies this effect. In this case, galaxies with rela- 
tively large optical luminosities but very low Lx (e.g. the recent 
detection by O' Sullivan & Ponman 2004) can be explained as 
objects still in the process of assembling the X-ray emitting 
ISM. Observations (e.g. O'Sullivan et al. 2001b; Samsom et al. 
2000) showing a correlation linking the age of the galaxies to 
their Lx, in the sense that younger systems exhibit lower lumi- 
nosities, give further support to this view. 

On the other hand, we cannot present more quantitative re- 
sults, because of the number of parameters involved. The as- 
sumed radius, in fact, can play an important role as well, be- 
cause the change of the dimension of the system leads to a 
change in the mean density and, thus, strongly affects the cool- 
ing (see Sec. 2.2.3). In Fig 6, the prediction for model Hal 
with R = 8/? e ff is plotted (empty hexagon). This galaxy ex- 
hibits a very high Lx, when compared with other models at 
nearly the same Lb, owing to its stronger cooling which allows 
more gas retention. When smaller radii (e.g. R = 3/? e ff) are as- 
sumed, however, the cooling becomes so strong that the galac- 
tic wind is delayed and lasts for a shorter time than in model 
Hal . Therefore, the models undergo a late time behaviour simi- 
lar to model HalO, namely an ISM gas mass of ~ 0. 1 -0.2Mi um , 
very low temperatures (~ 10 4 keV), and [Fe/H]> 1 .2, clearly at 
variance with observations. Since these models were run keep- 
ing all the other parameters fixed only for comparative pur- 
poses, they are not meant to match the optical properties and 
will not be discussed further. 

The natural conclusion is that the late accretion can play a 
role in the present-day L B - L x relation, its scatter being possi- 




Fig. 6. Prediction for the Lb - Lx relation for our model galax- 
ies compared to the best fit of O'Sullivan et al.'s (2001a, solid 
line) data, excluding bright cluster and group ellipticals as well 
as galaxies hosting an AGN. The dotted lines give a mea- 
sure of the scatter. Filled triangles: Lai, Hal. Open triangles: 
Lai, Hal in which we fix [Fe/H]=0. Stars: Models L and H 
with M acc = 0.007M lum . Filled squares: Models L and H with 
M acc = 0.005Mi um . Open squares: HO and L0. Hexagon: Hal 
with R=8fl eff . 

bly related to the differences in the gas reservoir from which a 
galaxy can form and possible different histories and efficiencies 
governing the secondary accretion. 

We anticipate from Sect. 5.4 that our model exhibits [Fe/H] 
higher than solar, at odds with observations. Therefore, in or- 
der to check whether a possible strong contribution from iron 
line emission can alter our results, we modelled two additional 
spectra by forcing [Fe/H]=0 in the ISM of both galaxies Lai 
and Hal. The luminosities obtained from these spectra (shown 
by empty triangles in Fig. 6) are lower than those emitted by 
Lai and Hal, and the differences are of the same order of those 
produced by a change in the accreted mass. 

In order to test the assumption of primordial composition 
for the second accretion flow, we run models Hal in which the 
accreted gas has Z = Z and Z = 0.3 Z , but always solar 
abundance ratios. From the entries in Table 3 it is clear that 
our results are robust against changes in the composition of the 
accreted gas. We did not comment cases with a Z = 0.3 - 1 • 
Z metallicity and a higher accretion rate, because the model 
predictions still show too cold galaxies to emitt in the X-Rays. 



5.3. The L x - T x relation 

In Fig. 7 we compare our predictions in the Lx - Tx plane with 
the data by Matsushita et al. (2000). In this particular case we 
show L x in the passband 0.5-10 keV, in order to be consistent 



Pipino et al.: X-ray properties of massive elliptical galaxies 



11 




-0.6 -0.4 -0.2 0.2 0.4 



Log kT x (keV) 

Fig. 7. Prediction for Lx - Tx relation for our model galax- 
ies Lai and Hal (filled triangles) compared to the data of 
Matsushita et al. (2000, squares with error bars). For compari- 
son we show the models Lai and Hal in which we fix [Fe/H]=0 
(open triangles). 




Fig. 8. Prediction for Lx— [Fe/H] x relation for our model galax- 
ies Lai and Hal compared to the data of Matsushita et al. 
(2000). The symbols are the same as in Fig. 7. Empty trian- 
gles refer to models Hal and Lai with the same Lx, but with 
[Fe/H]~ 0. 



with the observed energy range. The agreement is good for the 
Lai and Hal models, whereas the ones without accretion ex- 
hibit a slightly higher temperature. In any case, a factor of two 
in the temperature can be partly explained with variation in gas 
density and accretion history, but also might be related to the 
adopted one-zone formalism. Ferreras et al. (2002) presented a 
similar model for feedback and found a good agreement for the 
Lb - Lx relation, even though their low mass model exhibits a 
temperature of a few keV. They did not model any interaction 
with the environment, and thus any secondary accretion; more- 
over the fraction of gas allowed to escape was chosen a priori. 
This last point may represent another way to regulate the tem- 
perature in the ISM, without requiring secondary accretion. In 
fact, if we do not link the mass flow to the energetic budget and 
we have a situation in which AM/M gas < AE/E t h, the result is 
that the galactic wind removes the gas inhabiting the high en- 
ergy tail of the temperature distribution function, thus allowing 
the galaxy to retain more gas and eject more energy with re- 
spect to the case with no accretion. Therefore, the heating term 
cannot increase without limits even in presence of a continuous 
wind. 

We did not further explore this route, since it introduces 
an additional parameter to our formalism, but we stress that it 
might represent a natural solution in a more detailed modelling 
(e.g. two-phase ISM). 

A better agreement is achieved when we consider the mod- 
els Lai and Hal in which a [Fe/H]=0 ISM abundance ratio is 
forced (empty triangles in Fig. 7). 



5.4. The chemistry of the ISM 

Concerning the study of the chemical composition of the ISM, 
we first analyse the iron abundance, which is the best deter- 
mined among all the metals. Predictions for the L x - [Fe/H] x 
relation for our model galaxies Lai and Hal compared to 
Matsushita et al. (2000, squares) data are shown in Fig. 8 with 
filled triangles. As expected from the analysis of the fake spec- 
trum, our models predict a very high iron abundance, even if we 
compare our results with recent iron abundance determination, 
going from the typical ~ 1 - 2 Fe in the central region of giant 
elliptical galaxies (e.g. Xu et al. 2002; Buote 2002) up to the 
~ 3 - 4 Fe in the case of NGC507 (Kim & Fabbiano 2004). In 
the latter case, however, our model predictions are only within 
a factor of 2 from the observations. 

This high abundance is a clear consequence of our model, 
since all the iron produced, after the wind has stopped, is re- 
tained by the galaxies. In other words, the iron discrepancy 
seems to be still persisting. 

Among the possible physical mechanisms often invoked to 
solve this discrepancy (see Arimoto et al. 1997 for a compre- 
hensive analysis), we recall: (i) differential winds powered by 
SNe la and (ii) iron hidden in a colder phase. According to 
the former explanation, although our models stop the global 
wind phase, SNe la ejecta might be still able to escape the 
galactic potential well, thus lowering the iron abundance in the 
ISM (e.g. Recchi et al. 2001). In order to preserve the observed 
L x - L B and L x - T x relations, this could be obtained by e.g. re- 
quiring that AM ~ M 1SM holds at each time-step, namely only 
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the iron can escape, whereas the bulk of the gas is retained. 
This, in turn, implies that AM/M gas < AE/E^. 

The latter possibility, on the other hand, explains the low 
observed [Fe/H] by means of condensation in a colder phase 
(e.g. Fujita et al. 1996, 1997) or in dust grains. Recent obser- 
vations, in fact, claim that the dust mass in ellipticals could 
be ~ 10 6 ~ 7 M G , a factor of ten higher than previous estimates 
(Temi et al. 2004). If confirmed, this result will be at variance 
with the classical argument that the dust lifetime is short owing 
to the interaction with the hot ISM, and therefore it might rep- 
resent a viable solution to hide huge amounts of iron for a long 
period (Arimoto et al. 1997). 

However, recent observations (e.g. Bohringer et al. 1993; 
Finoguenov & Jones 2001; McNamara et al. 2000; Allen et 
al. 2001; Heintz et al. 2002) show that both the X-ray haloes 
and the region of the ICM very close to the central galaxy, are 
not uniform and that the presence of X-ray cavities, bubbles 
with relativistic gas that pushes away the 'thermal' gas, could 
be related to the AGN activity (Churazov et al. 2001). Central 
emission and the creation of radio lobes can drive a mixing be- 
tween the ICM and the ISM in such a way that part of the Fe 
is removed from the central galactic region. If we allow for the 
excess mass of iron to be replaced by primordial gas, our mod- 
els do not change their ISM gas mass, thus roughly preserving 
their L x , whereas the iron abundance is reduced to the observed 
values (empty triangles in Fig. 8). For example, in the simplis- 
tic case in which the metallicity is composed only of iron, if we 
assume that an AGN engine causes mixing to take place from 
the time at which the galactic wind stopped until the present 
day, and we want to remove the iron mass in excess (~ M' F s e M ) 
in order to end up with [Fe/H]= in the ISM, the rate of mass 
exchange is only ~ M^ e M ISGyr ~ 0.5MQyr~'. In the case of 
mixing involving non-primordial gas (e.g. [Fe/H]~ -0.5 in the 
ICM), instead, it should be either more extended in time or re- 
quire a slightly higher rate of mass exchange. Moreover, an 
additional source of heating (such as an AGN) might balance 
the strong cooling in model Ha 10, thus allowing the galaxy 
to maintain a longer wind phase. This, in turn, implies that 
more iron could be injected into the ICM and that the iron still 
present in the halo is diluted by a higher amount of gas with 
respect to our fiducial case (model Hal). At the same time the 
galaxy could exhibit the right temperature and gas mass in or- 
der to match the Lx - Lb relation. Moreover, the presence of 
an AGN could supply the amount of energy needed to mantain 
the galactic wind, in a scenario in which we onsider esraa = 0.1 
fixed for the whole galactic lifetime. In concluding, we note 
that evidences on the non-stellar origin of either a part of or 
all the ISM are given by the observation of subsolar iron abun- 
dance in galaxies which lie well below the average relation in 
the L x - L B plane (e.g. O' Sullivan & Ponman 2003), although 
the poor statistics might bias the observations (Matsushita et al. 
2000). Further aspects of this issue will be investigated in the 
future. 

Another possible solution can be a variation either with 
time or metallicity of the fraction of binary systems (e.g. De 
Donder & Vanbeveren 2002). This can be implemented in our 
code by assuming that the parameter A in the definition of SNe 
la rate (Matteucci & Greggio 1986) is not constant in time. 



We tested a model in which A scales linearly from to the 
assumed value of 0.09 with the gas metallicity, until the solar 
value (Z=0.02) is reached, and we did not find any significant 
reduction in the SNe la rate and the late time iron production. 
This is due to the very intense SF history experienced by ellip- 
ticals, which leads to a very fast increase of Z to the solar value 
in the first ~ lOOMyr of evolution. 

Concerning other chemical species, they represent a weaker 
constraint to our modelling than Fe, owing to the large observa- 
tional uncertainties still affecting their determinations. A sam- 
ple of recent observations is shown in Table 4. The galaxy iden- 
tity is given in Col. 1. Chemical abundances are presented in 
Cols. 2-9. In particular, subscript c means values measured in 
the core of the galaxy, whose radius is given in Col. 10, whereas 
subscript o is related to properties of the outskirts. Finally, Cols. 
11 and 12 show the adopted solar composition and the refer- 
ences. 

As it can be noticed from Table 4, there is a lot of confusion 
created by different adopted solar compositions. In the follow- 
ing we present our results adopting the solar meteoritic values 
of Anders & Grevesse (1989). Even though the presence of ra- 
dial gradients in the [a/Fe] abundance ratios is still debated 
(e.g. Gastaledello & Molendi 2002; Matsushita et al. 1997), 
the observations show subsolar ratios, especially in the cases 
of [O/Fe] and [Mg/Fe]. In our predictions, [Mg/Fe] traces very 
well [O/Fe], being both ~ -0.6 dex, whereas Si and Ca show 

a milder depletion ( 0.3 dex) with respect to the solar level. 

In any case, these predictions are ~ -0.2 dex lower than the 
measured values in the central (i.e. ~ 50kpc) regions of X-ray 
haloes, although still within their quoted ~ 1 - 2<x. 

On the other hand, the predictions from models undergo- 
ing a continuous wind (L0 see Table 3) are in better agreement 
with the observations presented in Table 4. It is worth noticing 
that these predictions reflect the high [a/Fe] ratios typical of the 
ejecta of low mass stars formed at the beginning of the galactic 
evolution and dying after 10 Gyr. This implies that, in order to 
reconcile our predictions with observations, a mechanism capa- 
ble of ejecting iron in a more efficient way with respect to other 
elements might be better than a simple dilution by primordial 
gas, which preserves the abundance ratios. 

We stress that the disagreement can be partly alleviated (at 
least at larger radii) by means of a multi-zone formulation, 
which takes into account the fact that at radii > \R e s the SF 
lasts for a much shorter timescale with respect to the galactic 
core, thus reducing the iron production and exhibiting much 
higher [a/Fe] ratios in the ejecta of low mass stars dying in the 
late stages of the galactic evolution. Therefore, more detailed 
modelling is required to assess this issue and to check whether 
the observed gradients can be reproduced. 

The observed ratios and the different behaviours shown by 
a-elements (and similar abundance patterns observed in the 
ICM, see next section) were used in recent works (Finoguenov 
et al. 2002, Gastaldello & Molendi 2002) in order to infer infor- 
mation on a possible change in the SNe II/SNe la ratio as well 
as on different explosion mechanisms as a function of galactic 
radius. Our opinion is that this kind of diagnostic is still prema- 
ture, given the uncertainties in the models used to fit the spec- 
tra (Mathews & Brighenti 2003; Gastaldello & Molendi 2002). 
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Galaxy 


[Fe/H] 


[O/Fe] 




[Mg/Fe] 


[S/Fe] 


[Si/Fe] 


[Ne/Fe] 


[C/Fe] 


[N/Fe] 


core 
(kpc) 


ab. 


ref. 


NGC5044 


-OJ-0.36,, 


-0.48 




-0.06 


-0.27 


-0.08 








48 


GS 




M87 


-0.11 


-0.44 




-0.15 






-0.26 





-0.21 




AG 


b 


M87 




-0.52,/ 


- 0.28,, 


-0.28 


-0.08 


C /0.1„ 


-0.52J-0.21,, 






2.5 


GS 




M87 


-0.009,/ -0.39,, 


-0.3V 


- 0.20,, 




0,/-0.09„ 


0.05,/ -0.04,, 










F 


f 


NGC1399 


~ 0.3,/ ~ -0.3„ 


















20 


AGm 


d 


NGC4636 


-0.06 


-0.24 




-0.13 






-0.09 




0.07 




AG 


e 



Solar abundance used: AG: Anders & Grevesse (1989, photospheric), AGm: Anders & Grevesse (1989, meteoritic), GS: Grevesse & Sauvall (1998), 
F: Feldman (1992). 

References: ": Buote et al. (2003), b : Sakelliou et al. (2002), c : Gastaldello & Molendi (2002), <*: Buote (2002), ': Xu et al. (2002), f\ Matsushita et al. (2003). 
Core: radius of the region quoted as 'galactic core' in the above-mentioned paper. 
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Fig. 9. Mass and energy flow rate versus time during the galac- 
tic wind for a Lai model galaxy. 

Moreover, we stress that it is incorrect to draw firm conclusions 
from observed abundance patterns just by relating abundance 
ratios to stellar yields and thus implicitly assuming the instan- 
taneous recycling approximation. 

6. The ICM 

Before concluding, we dedicate a section to the chemical en- 
richment of the ICM, since the improved chemical formalism 
presented in this paper allows us to follow the development of 
the galactic wind in a more detailed way with respect to pre- 
vious works (e.g. Matteucci & Gibson 1995; Martinelli et al. 
2000; P02). In Fig 9 we show, taking model Lai as an ex- 
ample, the temporal behaviour of both mass and energy flow 
rate. After a very short period in which they exhibit quite high 
values, most of the wind phase is governed by their slow de- 
crease as a consequence of the secondary gas accretion and 
the secular decline of SNe la. The mass flow rate has the 
same order of magnitude (in general within a factor of 2) of 
the value taken by the SFR just before the galactic wind, and 



this result is in agreement with the indications coming from 
the observations of starburst galaxies (e.g. Heckman 2002). 
On the other hand, the luminosity of the wind, after the very 
first intense blowout, is ~ 10 41 ergs _1 , giving a wind speed 
v ~ (2dE/dM) 1/2 ^ 200 km s _1 , again in agreement with the 
the typical estimates for starbursts (Heckman 2002). It has been 
claimed (Moretti et al. 2003) that, in the case of a Salpeter 
(1955) IMF, the SNe la-driven winds can last for a long pe- 
riod, in order to inject into the ICM a suitable amount of iron, 
only in the case of very high SN thermalization efficiency (i.e. 
esNia = 1, P02). This is not evident, since in this paper we show 
that the wind phase is about 1/3 of the galactic lifetime, being 
fully supported by SNe la activity, even if esNia *£■ esNia,P02 = 1 
for most of the galactic evolution. 

Assuming that the metals observed in the ICM come from 
elliptical (and perhaps SO) galaxies (Arnaud et al. 1992) and 
following the same procedure as described by P02, we found 
no significant changes in the predicted amounts of gas, iron, 
oxygen and energy ejected into the ICM by the model galaxies, 
with respect to the old model. In particular, the total amount of 
oxygen ejected into the ICM, 0.43 x 1O 9 M (0.44 x 1O 1O M ) by 
model Lai (Hal), is very similar to the entries of P02 Table 1, 
whereas the new models lose 0.28 x 1O 9 M (0.23 x 1O 1O M ) of 
iron, which is a factor of 2 higher with respect to P02 model I 
(but in agreement with their estimate for model II), owing to the 
more extended SF history which has increased the amount of 
iron released promptly at the beginning of the wind. Therefore, 
even though the new models predict ellipticals slightly younger 
than those in previous works, the amount of iron in the ICM 
can still be explained by means of the galactic winds, and we 
can recover the typical ~ 0.3 Fe & (e.g. recent review by Renzini 
2003; White 2000) by assuming that most of the gas in the ICM 
has a primordial origin (Matteucci & Vettolani 1988, David et 
al. 1991, Renzini et al. 1993, P02). Moreover, the bulk of iron 
is ejected before at z > 1, in agreement with recent observations 
on high redshift galaxy clusters (Tozzi et al. 2003). The amount 
of energy ejected by our model galaxies into the ICM, is simply 
the integral of the curve in the lower panel of Fig. 9. The results 
are very close to those derived by P02 for their Model I, there- 
fore still not sufficient to provide the required ~ IkeV in order 
to break the cluster self-similarity (e.g. Tozzi & Norman 2001; 
Loewenstein 2001; Bialek et al. 2001; Borgani et al. 2001). 

Concerning the oxygen and other a-elements, new data 
appeared in the last few years (see the recent review by 
Loewenstein 2003) and a more consistent picture seems to be 
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emerging from those measurements. In particular, more accu- 
rate [O/Fe] ratios are observed for a larger sample of galaxy 
cluster (Peterson et al. 2003) and, added to previous detections 
(e.g. Tamura et al. 2001), seem to point toward generally sub- 
solar values. The typical values lie around [O/Fe] 0.2 (e.g. 

Peterson et al. 2003), but the dispersion is quite high, so that 
our predictions ([O/Fe] 0.6) are within ~ 2cr (given the typ- 
ical uncertainty quoted by Peterson et al. 2003). Therefore, we 
stress that the new observations allowed by XMM and Chandra 
satellites seem to converge toward a better agreement with our 
predictions. On the other hand, other a-elements exhibit differ- 
ent degrees of enhancement with respect to iron (Baumgartner 
et al. 2004), although the uncertainties associated with their 
observations are still very high. In particular, our predictions 
for Ca ([Ca/Fe]~ -0.4) agree with the values measured by 
Baumgartner et al. (2004), whereas we find [Mg/Fe]~ -0.7, 
which is at variance with the typical [Mg/Fe]~ (e.g. Ishimaru 
& Arimoto 1997; Peterson et al. 2003, but note that their Icr 
error is ~0.4 dex). In any case a larger sample of clusters is re- 
quired before drawing strong conclusions also on a-elements. 
Furthermore, a new kind of theoretical modelling should be de- 
veloped, since in many cases the abundance ratios as well as 
the temperatures of the X-ray haloes of bright central galax- 
ies seem to vary with radius until they reach the typical values 
of the ICM. For example Peterson et al. (2003) explicitly no- 
ticed that their results resemble those of Xu et al. (2002) for 
NGC4636. At the same time, the galaxies located at the center 
of galaxy clusters, show properties that seem to correlate with 
the excess of iron measured in the center of cold core clusters 
and the ICM temperature (De Grandi et al. 2004). 

Finally, we stress that if the mechanisms invoked to solve 
the Fe discrepancy in the ISM were at work, our conclusions 
on the ICM would not have changed. For istance, the Fe ex- 
cess in the ISM which should be removed from the halo in 
order to achieve a solar abundance in the X-ray spectrum is 
0.29 x 1O 9 M for model Hal. Therefore the mass of Fe ex- 
pelled in this way is only the 10% of the total mass of Fe al- 
ready ejected into the ICM during the galactic wind, and also 
the effects on the predicted ICM abundance ratios will be neg- 
ligible. Among the suggested mechanisms the presence of dust 
could be relevant and since different scenarios could be at work 
together, we believe that the real effect of Fe removal from the 
ISM to the ICM, will be even lower. 



7. Conclusions 

The model presented in this paper is a first step in the self- 
consistent study of both optical and X-ray properties of ellipti- 
cal galaxies by means of a chemical evolution code. In order to 
do that we updated previous chemical evolution codes (PM04) 
by implementing a new energetic treatment which takes into 
account for the first time cooling and heating processes occur- 
ring in the ISM. Adopting the procedure described by Kawata 
& Gibson (2003b) we are able to predict the X-ray spectrum 
of our model galaxies taking into account the presence of line 
emission and the actual energy range used in the observations. 
Here we present a summary of our main conclusions: 



- Our new models confirm the previous finding of PM04, 
namely that SF and infall timescales decreasing with galac- 
tic mass are needed to explain the optical properties of el- 
liptical galaxies. 

- At the same time we reproduce the L x - L B relation in the 
ISM of bright ellipticals. Our model can also predict the 
right ISM temperature in cases in which a mild accretion 
from the surroundings occurs, and within a factor of 2 from 
the observed temperature at a given L x for models with- 
out accretion. The synthesised X-ray spectrum is in a good 
qualitative agreement with the typical observed ones. 

- In order to model in a more realistic way the formation 
history of ellipticals, we implemented the possibility of a 
secondary accretion episode spread over the entire galactic 
lifetime. This mild inflow seems to play a non-negligible 
role in governing the late-time behaviour of the ISM, in par- 
ticular in driving L x - L B and L x - T x relations, whereas 
its effects on the colours and the line-strength indices are 
negligible. 

- Variations in the ISM density (adopted radius) and the 
amount of mass accreted during the secondary infall might 
account for the large spread in the observed L x - L B rela- 
tion. 

- Our model predicts [a/Fe] abundance ratios in the ISM 
which are in agreement with the data coming from very 
recent observations with XMM and Chandra satellites. 

- The iron discrepancy between model predictions and ob- 
servations still persists, probably indicating the presence of 
other mechanisms acting at galactic scales. We suggest that 
the mixing of gas driven by AGNs can preserve the gas 
mass (and thus the X-ray luminosity) while diluting the iron 
abundance. 

- The new energy formalism implemented in the chemical 
evolution code allows us to follow in a more detailed way 
the evolution of mass and energy flow into the ICM with 
respect to previous work (e.g. Matteucci & Gibson 1995, 
Martinelli et al. 1998, P02). Nevertheless the predicted 
amount of Fe, [a/Fe] ratios as well as the energy injected 
into the ICM are very similar to P02 results. Therefore, we 
confirm that SNe la are fundamental in providing energy 
and iron to the ICM and we notice that very recent obser- 
vations have reduced the gap between our predicted subso- 
lar [a/Fe] and the measured ones. Even though our fiducial 
case has esNia <T the SNe la activity can power a galactic 
wind lasting for a considerable fraction of the galactic life- 
time and reproduce the observed iron abundance in clusters. 

It is important to stress that we achieved these results in the 
assumption of a one-zone formalism with a single gas phase in 
which the model galaxy is treated like an object with uniform 
density. Future developments of this kind of models should in- 
clude the gas density profile, the possible presence of an AGN 
and a multi-phase approach. In this way we will obtain a more 
reliable estimate of the ISM temperature and X-ray luminosity, 
which are strongly related to the gas profile, and we might be 
able to model the gradient in the abundance ratios observed in 
both the ISM and the region of the ICM close to the brightest 
cluster galaxies. 
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